Disentangling Strong Dynamics through Quantum Interferometry 
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We present a new probe of strongly coupled electroweak symmetry breaking at the 14 TeV LHC 
by measuring a phase shift in the event distribution of the decay azimuthal angles in massive 
gauge boson scattering. One generically expects a large phase shift in the longitudinal gauge boson 
scattering amplitude due to the presence of broad resonances. This phase shift is observable as an 
interference effect between the strongly interacting longitudinal modes and the transverse modes of 
the gauge bosons. We find that even very broad resonances of masses up to 900 GeV can be probed 
at 30 significance with a 3000 fb~* run of the LHC by using this technique. We also present the 
estimated reach for a future 50 TeV proton-proton collider. 


Introduction. One of the most important goals of the 
Large Hadron Collider (LHC) is to find the nature of the 
mechanism of electroweak symmetry breaking (EWSB). 
The discovery of the 125 GeV Higgs-like object [I [2] is 
a major milestone in this direction. However, in order 
to truly understand EWSB, we would like to learn the 
origin of the longitudinal components of massive elec- 
troweak gauge bosons (V). Broadly speaking, models of 
EWSB fall into two categories, those where the longitudi- 
nal components of the Vs are a) weakly interacting or b) 
strongly interacting. The most popular examples in the 
first category are the Standard Model (SM), with one or 
more elementary Higgs multiplets [3], and its supersym- 
metric extensions [4]. In the second category, one or more 
strongly interacting sectors appear at the TeV scale and 
are responsible for EWSB. The Higgs-like boson and the 
longitudinal Vs could arise as pseudo-Nambu-Goldstone 
bosons (PNGBs) in this scenario [2I]. A definite dis- 
tinction between these two cases at the LHC would be 
very important and serve as one of the first crucial steps 
towards a full understanding of EWSB mechanism. 

One universal consequence of a strongly coupled 
EWSB sector is an enhancement of the longitudinal 
gauge boson scattering amplitude at high energies over 
the standard model expectation [22]. In order to measure 
such an effect, one would look for an enhanced rate for 
total gauge boson scattering at high invariant mass for 
the gauge boson pair, and then use the polar angle dis- 
tribution of the gauge boson decay products to measure 
their longitudinal polarization fraction [23]. 

In addition to an enhancement of the magnitude of 
the amplitude, one also expects a large phase shift in the 
amplitude, as in the case of scattering through a reso- 
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nance. In this paper, we would like to seek strategies to 
experimentally probe this phase shift induced as a con- 
sequence of strong dynamics. We will use the azimuthal 
angle correlations of the Vs’ decay products and show 
that the strong phase shift shows up as a modification 
to the interference effect between the longitudinal and 
transverse V polarizations in the azimuthal angle dis- 
tributions. Thus, the phase shift is turned into an ob- 
servable that can be used as a complementary probe (in 
addition to the longitudinal V scattering rate) of EWSB 
from strong dynamics. 

PNGB scattering and parameterization of the 
phase shift. Pion scattering provides a realistic ex- 
ample of strongly coupled PNGB scattering. By look- 
ing at the experimental 77 scattering data [9], we can 
see that a large phase shift exists in the form-factor of 
both mz — am and ete~ — nr. The phase shift ô in 
ete” — nr scattering versus energy y's is shown by the 
black data points in Fig. {1} We can see that the ampli- 
tude undergoes a large phase shift when vys is near the 
mass of the p meson (760 MeV). 

The low energy effective theory of pion scattering is 
described by a chiral Lagrangian. In this description, 
the scattering amplitude in the isospin (J) and partial 
wave (J) channel of nr — rr is given by, 


My x sin(577)e?"/ = cry el te, (1) 
T 


where fr = 84 MeV is a low-energy constant of chiral 
perturbation theory (which differs from the physical pion 
decay constant, fPPYS = 92 MeV that arises away from 
the chiral limit), and the last equality follows from the 
low energy theorem prediction of the behavior of the scat- 
tering amplitude with s, which also determines the con- 
stants cry. Naive extrapolation of this form of the am- 
plitude violates unitarity for energies ,/s comparable to 
the cutoff 47 f, of the effective theory. In particular the 
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FIG. 1: The phase shift 5 in ete” — mr scattering versus 
energy ys. The experimental data (black points) are from 
Ref. [9]. The red solid line is the ansatz in Eq. we have 
used to fit the phase shift using the parameters m = 760 MeV 
and I = 140 MeV. Using dispersion relations we can relate 
the phase shift to the magnitude of the form factor |F(s)|. 
For P(s) = 1, the magnitude is shown by the blue dashed 
line. 


low energy theory fails to predict a resonant enhance- 
ment of the magnitude of the amplitude, as well as the 
large phase shift expected from the exchange of (J = 1) 
vector-meson resonances. 

There are several ways to incorporate the effect of the 
vector resonances into the scattering amplitude calcula- 
tions of the low energy theory. One approach is to add 
in a broad vector resonance by hand to the theory. How- 
ever, we will choose to adopt a different approach that 
manifestly maintains the unitarity of the amplitude at 
very high energies and emphasizes the central role of the 
phase shift. 

We will choose to multiply the tree level amplitude (in 
the J = 1 channel) by a complex form factor F(s). The 
entire form factor can then be extracted from its phase by 
using analyticity arguments to define an Omnès function 
[8] and assuming no inelastic channels. For a given form 
factor, F(s), applying the subtracted dispersion relation 
to log(F'(s))/s, we have 

co 

F(s) = P(sjexr(— | asse = - 5 
T Jo s'—s—te s 
When the phases 6(s’) goes beyond 27 (for instance, mul- 
tiple resonances with additional branches), the additional 
27n phase factors can be recast into P(s) as a polynomial 

factor with P(0) = 1. 

For simplicity, we use an ansatz for the phase, 


5(s) = ArcTan[sl’'/m(m? +T? — s)], s < m? 3 
(s) = ArcTan[s’'/m(m? +r? — s)] +7, s >m?, (3) 


which approaches a constant at high energy according to 
unitarity. We can see that this ansatz can fit the phase 
of the mr scattering data very well, as shown in Fig. 


Nevertheless, our parametrization is general and the does 
not rely on the specific form of strong dynamics. From 
Eq. (2), we can construct the full form factor from its 
phase 6(s), 
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F (s) = P(s) (4) 
We can see the behavior of |F (s)| in Fig. |1| At small s we 
are far from the resonance and the form factor is unity. 
However as we approach the resonance the magnitude 
of the form factor grows large. For large values of \/s 
beyond the resonance, the form factor falls off rapidly 
with energy. 

Longitudinal weak bosons scattering. In a strongly 
interacting EWSB theory the longitudinal components 
of the weak bosons can be approximately regarded as 
PNGBs, and their interactions can be described at low 
energies by a chiral Lagrangian. Thus, naive extrapola- 
tion of the scattering amplitude would lead to a similar 
problem of unitarity violation. We can expect a similar 
resolution to this problem as in the case of pions, where 
new resonances unitarize the scattering amplitude. 

However, there is one key difference between weak bo- 
son scattering and pion scattering. Namely, we have al- 
ready discovered a 125 GeV (Higgs-like) scalar object 
that couples to weak bosons. The exchange of this scalar 
object will partially restore unitarity in the longitudinal 
V scattering amplitude. If the couplings of this scalar 
object to Vs is exactly the Standard Model value, corre- 
sponding to a non-composite Higgs boson, then we would 
have complete unitarity restoration just from exchange of 
this scalar object. 

For composite Higgs models, we can use the general- 
ized Adler-Weinberg sum rule (in the limit of vanishing 
gauge couplings) [10] to relate the Higgs coupling to Vs 
to an integral sum of longitudinal gauge boson scattering 
cross-sections in various isospin channels, 


co 1.2 

1-a = | ZE (20s) + 30 (s) — 501,18). (8) 
0 8 

Here, a parameterizes the ratio of the Higgs boson cou- 

pling to Vs over its SM value. From the latest fits [II] 

we can see that 0.8 < a < 1.2 at the 2o level. 

For a not equal to 1, we need additional contributions 
from strong dynamics to restore unitarity in longitudinal 
gauge boson scattering. For simplicity, we will also as- 
sume vector meson dominance for the remaining partial 
unitarity restoration, as in the case of pion scattering, 
restricting the new physics contribution to the J = 1 
channel. 

We can now parameterize strong dynamics in longi- 
tudinal V scattering by introducing a form factor that 
multiplies the amplitude just as we did in the case of a7 
scattering. We will assume the same ansatz for the phase, 
since it maintains unitarity manifestly. The role of the 
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FIG. 2: The currently excluded region of form-factor parame- 
ters (shaded region) using constraints from CMS searches for 
W' + WZ resonances [12]. (The ATLAS bound is slightly 
weaker and is not shown.) 


composite Higgs boson in unitarization can be accounted 
for without explicitly including the Higgs exchange dia- 
grams. Instead, we note that J = 1, longitudinal V 
cross-section must be rescaled by a factor of 1 — a? (com- 
pared to the case where no Higgs-like boson is present). 
Thus, we can simply rescale our form factor from Eq. 
by v1 — a? to account for the presence of the 125 GeV 
boson. 

Assuming a = 0.8, we show the currently excluded 

region for different form-factor parameters m and I in 
Fig. |2| using the latest LHC searches for W’ => WZ 
resonances [12]. This search strategy only probes the 
enhancement to the total rate for WZ production, or 
equivalently it is a probe of |F (s)|, but it does not probe 
the phase shift directly. In addition, the search uses a 
narrow WZ invariant mass window to suppress back- 
grounds and is therefore insensitive to very broad res- 
onances T/m > 20%, that are characteristic of strong 
dynamics. 
Observation of Longitudinal V scattering. There 
are two methods to observe VV scattering at the LHC. 
One is the widely used weak boson fusion process pp —> 
VV jj, where the two forward jets can be used to sup- 
press the large SM backgrounds. The other method is 
the rescattering process pp + VV, which is not consid- 
ered to be a promising channel for most studies due to the 
large SM backgrounds. However, this channel is not sup- 
pressed by the small effective-V luminosity and it also 
has better access to higher energies of the VV system. 
Since we are trying to observe the azimuthal angle cor- 
relation which arises from a quantum interference term, 
we do not have to suppress the SM backgrounds and we 
will consider the rescattering process in this paper. 


Angular correlation. Let us consider W* Z production 
from a ud initial state. We will consider a modification of 
the SM amplitude of longitudinal W* Z production in the 
J = 1 channel by a complex form-factor V1 — a? F(s), 


FIG. 3: The kinematics of Wt Z production and decays. 


where F(s) is given by the form in Eq. (4). We will 
study this process at high energies where both the W+ 
and the Z decay leptonically. The kinematic dependence 
of the production and decay amplitudes are as follows: 
My: u(k1, —) d(k2, +) => W* (qi, A1) Z (qo, 2); Mo : 
W*(m,A1) = v(p1, —) I+ (pe, +) and Ms : Z (qa, A2) = 
l3 (p3, h) I3 (pa, —h). The parameters in the parentheses 
are the particle momentum and helicity respectively. The 
phase space of this process has five independent angles 
which are defined in the center-of-momentum frame as 
follows: the production angle (©), two polar decay angles 
(01, 82) and two azimuthal decay angles (%1, ¢2) which 
can be though of as the rotations of the W/Z decay 
planes (Aw ,fz) about the Wt/Z momentum axis, and 
are measured relative to the production plane ù. The 
three planes are defined as ù ~ kı X qı, NW ~ qı X p2 
and fiz ~ p4 X q2. All these kinematic variables are 
presented in Fig. 

The phase shift from the strong dynamics only af- 
fects the longitudinal-longitudinal combination of W+ Z 
modes and shifts the corresponding amplitude by an en- 
ergy dependent phase, ô. This phase shift will enter into 
the azimuthal angle correlation of the Wt/Z decay as 
an interference effect between the various polarizations of 
the vector bosons. To see this, recall that W*/Z decay 
produces an azimuthal angular dependence exp(is,@) in 
the amplitude, where s, and ¢ are respectively, the spin 
projection and the azimuthal angle rotation about the 
W* /Z direction. For a given interference term between a 
general helicity combination (A1, A2) and the (0,0) combi- 
nation, we find that the relevant terms in the differential 
cross-section are of the form |a,e"(191~ 202) + azet’ |? ~ 
cos(A11 — dodo + ô) D sin(A1¢1 = A2db2) sin ô [14]. 
Thus, the sin(A1¢; — Ag¢2) modes in the azimuthal an- 
gle correlation strongly suggest the existence of strong 
dynamics. 


The production amplitudes can be separated by the 
different W*Z helicity combinations (A;,A2). Among 
the nine different helicity combinations, there are three 
leading contributions from (—,+), (0,0) and (4, -). 
The four other significant ones are (—,0) ~ (0,+) and 


(+,0) ~ (0,—), all of which have a relative suppression 
~ mw /,/s compared to the leading modes. The (+, +) 
and (—,—) combinations are too small to affect the kine- 
matic distributions. We note that (a): (—,+)/(+,—) 
dominate as © — 0 or m because of the t-channel pro- 
duction. (b): Numerically the difference between (—, 0) 
and (0,+) or (+,0) and (0,—) is negligible. We param- 
eterize the production amplitudes as Mı(—, +) = A, 
M,(0,0) = Be, Mi(+,-) = C, Mi(-,0) = D, 
M,(+4, 0) = E, Mı (0, +) = F and M,(0,—) = G. 
Here, all amplitudes depend on the center-of-mass en- 
ergy and on O. Similar behavior has been pointed out in 
ete —> WtW [I5]. 

The full differential cross-section can be obtained from 


2 
D Lae My (9) M3 (1, 1) M3? (h, 02, 2) ’ 

where the W+, Z decay amplitudes are MÀ (01, ġ1) 
= gwlialdy, (0; )etA1?1, M3? (h, 01, Qı) = g'\qa\dy. (02 Te 
(h — 1)r/2)e—**242, Here gz ~ —gf are the couplings 
of the Z to different lepton helicities. The polar angle 
dependent function d4 (0) = Ja +cos 0), do(@) = sin 8. 
If we integrate over both polar angles 6; 2, there is an 
approximate cancellation in the sin(¢@, + ¢2) and sin(¢2) 
correlation between coso > 0 and cos62 < 0 due to 
g% © —gz. Therefore, we only integrate the differential 
cross-section over either cos #2 > 0 or < 0 to obtain: 


on 5 a 5((# + p A +G? r») + Bsind (3 (+4(A + C) — 3re( A — C)) sin(d1 + Q2) 
a dad HEO aO sin) +] . (6) 
| 
where e = ((97)? — (94)°)/((97)? + (g4 )°) ~ 0.22 and Mı 


the ellipsis refer to the interference terms with cos ¢-type 
dependence. The o stands for o(cos62 Z 0) and H is 
overall background H? = (A? + B? + C? + D? + E? + 
F? + G?). 

In Fig. [4}we plot the relative coefficients of the different 
sin ¢ correlations as a function of cos © for ys = 1 TeV 
using the expressions in Eq. (6). Measuring a non-zero 
coefficient for any of the sin ¢ modes is a positive indica- 
tor of strong dynamics. Integrating over the entire cos O 
range would lead to cancellations that would dilute the 
significance of the probe. Thus, an optimal strategy is 
to use a maximum likelihood analysis on the measured 
do+/dcos Odd, dd¢2 distribution in the data to look for 
all non-zero sin(A;¢1 — A2¢2) coefficients. 


Procedure. At the LHC, there are various kinematic 
ambiguities that must be incorporated into the analysis. 
We will choose to study only the sin ¢; mode, which will 
yield a much more transparent analysis that is robust to 
the kinematic ambiguities. 

We simulate the process pp > Wt Z — Itvltl- with 
the form factor v1 —a?F(s) from Eq. multiplying 
the J = 1 channel of the (0,0) helicity amplitude. We 
have taken a = 0.8 as a benchmark, which assumes that 
such a deviation in Higgs couplings would continue to 
be allowed with future data from the LHC. We then scan 
over different form factor parameters m and I, which will 
give rise to different phase shifts. Our event simulation is 
purely at the parton level which is sufficient for the fully 
leptonic final state. We used HELAS [16] to calculate 
the helicity amplitudes for the full process. LHApdf [I7] 


FIG. 4: The coefficients of the production amplitudes in front 
of the azimuthal sin ¢ correlations as a function of scattering 
angle cos@ at ys = 1 TeV from Eq. (6p. The solid blue, 
green, red and dashed blue, red lines stand for the coefficients 
in front of sin(¢?1 + ¢2) (cos 62 > 0), sin ¢1, sin 2 (cos 62 > 0), 
sin(d1 + ¢2) (cos 62 < 0), sin dz (cos 02 < 0) respectively. 


was used to fold in the parton distribution functions for 
the protons using the pdf set CTEQ6L [I8]. An adaptive 
Monte-Carlo package, BASES [19], was used to perform 
the integration over phase space and to study differential 
cross-sections. Our simulation approach to calculation of 
helicity amplitudes allowed us to insert the form factor 
specifically into the longitudinal gauge-boson scattering 
channel. 

The cuts used are: 1): Ar > 0.4 separation between 
leptons; 2): pr > 20 GeV and |n| < 2.8 cuts on the 


leptons. 3): Z-reconstruction cut: We require that two 
opposite sign leptons reconstruct to give the Z mass. 4): 
Missing Er cut > 20 GeV. 5): Invariant mass cut of the 
W/Z system between m +T [24]. 

When we reconstruct the events, there are two misiden- 
tification issues that arise that lead to a fourfold ambi- 
guity in the kinematics: a): The u-quark direction is 
unknown. b): There is a two fold ambiguity in the neu- 
trino momentum along the beam axis. First, consider a 
misidentification of the u-quark direction. This leads to 
misidentifying 0 > 17-0, ¢, > 7+41, $2 > T+. The 
azimuthal angle correlation, sin @,, is odd under such a 
misidentification. Note that from the solid green curve in 
Fig. [4] the coefficient of sin ¢; is also approximately odd 
under © — 7— 0. Thus, if we study the sin ġı mode 
for either 0.1 < cosO < 0.9 or —0.9 < cosO < —0.1 we 
find that it is robust to misidentifications of the u-quark 
direction. 

The presence of a false solution for the neutrino mo- 
mentum would distort the azimuthal angle correlations 
that we seek for ġı. However, to study the sin ¢; mode, 
we can simply measure the up-down asymmetry with re- 
spect to the ¢; = 0 (production) plane to find the sizes 
of such correlations. We will demonstrate that the up- 
down asymmetry is the same for both the true and the 
false solutions. 

For a given ¢; azimuthal angle correlation we have, 


do 


doy cos 820 


~ Ao + A, cos ġı + Ap cos 2¢; + Bı sind, . 
(7) 


We define the events going “above” the plane for sin ¢1 > 
0 and going “below” the plane for sind; < 0. Therefore, 
the up-down asymmetry can be defined as 


Ni — N 2 Bı 
AS| os 020 = Ni +N ~ at Ao ’ (8) 


where N/N_ are the number of up/down events respec- 
tively. 

The up or down events for ¢, can be defined by the 
sign of the scalar triple product 


SGN = sgn (î - p2) = sgn ((kı x q1)-p2) - (9) 


For a particular event if SGN > 0(< 0) then we increment 
N+(N-). The normal vector to the production plane 
A = ky X q = kı x (pı + p2) is independent of the v 
momentum along the u-quark direction and hence SGN 
is insensitive to the difference between the true and false 
solution. In addition to this, the asymmetry variable has 
the advantage of being insensitive to a number of cuts 
such as rapidity and pr cuts that would otherwise distort 
the angular distribution. 


Results. If the background fluctuation is Gaussian, the 
statistical significance of the nonzero asymmetry is given 


14 TeV LHC 50 TeV pp collider 
| T/m > |0.3[0.4][ T/m > [o.2[0.3[0.4 
m (GeV) m (GeV) 


800 [3.44.8 1000 |5.7/6.7|9.1 
900 [2.4|3.3 1100 |3.7|4.8|6.8 
1000 |1.7|2.4 1200 |2.4}3.4/5.1 


TABLE I: Table showing the cross-sections (in fb) for the 
fully-leptonic process pp > W*Z —> I*vl™l for different 
form-factor parameters in a composite Higgs model (assuming 
a = 0.8) and with the cuts described in the text (including 
the mass window selection). Cross-sections are shown for the 
14 TeV LHC and at a future 50 TeV pp collider. The typical 
cross-sections are of the order of a few fb. 


by, 


js N = |AS|VN, (10) 


where N = N+ +N- is the total number of events. As a 
rule of thumb, we found that choosing an invariant mass 
window between m +T seemed to optimize the trade- 
off between picking up a large |AS| by being close to the 
resonance, while still keeping a sizeable number of events. 

In Tab. |I| we show the cross-section for the process 
under consideration at the 14 TeV LHC and at a future 
50 TeV pp collider for different choices of form factors by 
varying over the parameters m and I. 

The form factors that we consider lead to typical asym- 
metries of the order of 5 — 10%. In Tab. we show 
the significance of the asymmetry measurement at the 
14 TeV LHC with 3000 fb~+ of data for different choices 
of form factors by varying over the parameters m and rF. 
We show the results using expected statistics including 
both the W*Z and W-Z fully-leptonic modes. 

We find that new wide resonances can be probed at 
the 30 level for masses up to more than 900 GeV. This 
further motivates an extended run of the LHC should an 
excess in W*Z be discovered. A future 50 TeV proton- 
proton collider could probe resonances up to around 1.2 
TeV with the same luminosity. An even higher energy at 
a future collider would make the sin(¢1 + ¢2) interference 
term the dominant piece and would require a different 
analysis strategy. 

There are several theoretical and analysis issues that 
could potentially increase these significances in a more 
sophisticated search and motivates an elaboration of our 
work. 1): In the sin ¢, mode search, it is possible to open 
up the hadronic decay modes of the Z with boosted tag- 
ging techniques [20]. 2): In addition, using a multivariate 
analysis and incorporating the sin@2 and sin(¢, + ¢2) 
modes could also bolster this result. 3): Allowing for 
I = 2 resonances implies that the vector resonance form 
factor could be scaled by a factor larger than V1 — a?. 
4): Multiple resonances in a narrow mass window could 


14 TeV LHC 50 TeV pp collider 
T/m > [0.3]0.4|| T/m > |0.2]0.3]0.4 
m (GeV) m (GeV) 

800 4.4|4.4 1000 |4.8|5.0|4.9 

900 [3.2/3.2 1100 13.313.6|3.6 

1000 |2.3|2.3 1200 |2.4|2.6|2.7 


TABLE II: Table showing the significance of the asymmetry 
variable, from Eq. {10}, using both the W~Z fully-leptonic 
modes for different form-factor parameters in a composite 
Higgs model (assuming a = 0.8). The results are shown for 
an integrated luminosity of 3000 fb~! at the 14 TeV LHC and 
at a future 50 TeV pp collider. 


also yield an enhancement in the longitudinal scattering 
cross-section which would show up as the P(s) factor 
mentioned earlier. 


Conclusions. We have proposed a novel technique to 
disentangle the dynamics of a strongly coupled EWSB 
sector by measuring a phase shift in the the decay az- 
imuthal angle correlations in massive gauge boson scat- 
tering. Our results show that a simple up-down asym- 
metry in leptons from W decay in pp > W~Z is robust 
to a number of event reconstruction ambiguities and is 
a good probe of broad resonances from strong dynamics. 
This strongly motivates a high luminosity run of the 14 
TeV LHC. A future 50 TeV pp collider could yield conclu- 
sive evidence of resonant behavior in the presence of an 
excess of WZ events at the LHC. Furthermore, we have 
outlined several analysis strategies and theoretical issues 
which would significantly increase the reach of searches 
based upon this technique and could lead to a promising 
signal at the next run of the LHC. 
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There could be other non-universal consequences, such 
as resonances, enhanced di-higgs production rate [6], etc. 
See Ref. and references therein. 

In practice this requires looking for a broad excess and 
placing an invariant mass cut on the W Z. We found that 
our results were not significantly affected by choosing a 
different invariant mass cut window, indicating that a 
precise optimization of this cut may not be very impor- 
tant. 


